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1-H Abstract 

in 

(^ The effect of edges and apertures on the Casimir energy of an arrangement of plates and 

■^ boundaries can be calculated in terms of an effective nonlocal lower-dimensional field theory 

<^ that lives on the boundary. This formalism has been developed in a number of previous papers 

• and applied to specific examples with Dirichlet boundary conditions. Here we generalize the 

1^ formalism to arbitrary boundary conditions. As a specific example, the geometry of a flat 

.^H plate and a half-plate placed parallel to it is considered for a number of different boundary 

/\ conditions and the area-dependent and edge dependent contributions to the Casimir energy 

^ are evaluated. While our results agree with known results for those special cases (such as the 

Dirichlet and Neumann limits) for which other methods of calculation have been used, our 

formalism is suitable for general boundary conditions, especially for the diffractive effects. 



1 Introduction 

The Casimir effect [1], the classic example of the influence of boundary conditions in quantum 
field theory, has been of considerable interest over the last several years. A vast body of 
literature has emerged with a number of new geometries being explored, both analytically and 
numerically and both at zero and nonzero temperatures p]. Recent developments in nano- 
machinery have also provided further impetus to these efforts. Geometries with edges and 
apertures are a particularly interesting set because of the possible interplay of diffractive effects 
and boundary conditions on the fields. The analytic calculation of the propagators in the given 
geometry with the given boundary conditions is rarely feasible and numerical and analytical 
approximations have to be used; diffractive effects are then a formidable computational task. 
We have recently developed a formalism which focuses on the diffractive effects [3]- [6]. In our 
approach, an effective nonlocal lower dimensional field theory defined on the boundaries, with 
a boundary action Sb, is first extracted. In the subsequent analysis of this lower dimensional 
theory, apertures and edges can be easily incorporated leading to a systematic calculational 
scheme for diffractive effects. The results, both for zero and nonzero temperatures, were in 
good agreement with alternate methods of calculation, such as the numerical world- line method 
[71 [8] , and the scattering matrix method P, HIOl] and in those geometries where such calculations 
had been done. We also obtained a universal low energy theorem for the Casimir force between 
holes on a plate [4j. 

In this paper, we extend this formalism addressing two important and related issues. The 
boundary action Sb is a functional of the field on the boundary and hence it is naturally 
defined for Dirichlet boundary conditions for the propagator or the fields. So the first question 
is: How do we apply our formalism to the case of Neumann or even more general boundary 
conditions? Although the case of two full parallel plates with no edge or diffractive effects 
for a scalar field with Robin boundary conditions has been studied before in |1H I12j . Our 
focus here is on the edge and diffractive effects. While Dirichlet, Neumann or Robin are the 
easiest to work with, the general boundary conditions, say on a Laplace-type kinetic energy 
operator, are given by the von Neumann theory of self-adjoint extensions characterized by a 
self-adjoint operator on the space of boundary values [13^ I14j . However, there is a problem 
with the simple and straightforward generalization to arbitrary boundary conditions because, 
generically, the (minus) Laplace operator has negative eigenvalues, suggesting an instability for 
the quantum theory |14j. On the other hand, physically, the boundary of interest is produced 
by the insertion of a material plate or some such object, and it would be rather strange if 
this process leads to a global instability. So the next question would be: How do we square 
these concerns? Further, it is not clear that, in a full hedged interacting quantum field theory, 
mimicking the material plate by boundary conditions is an adequate characterization. These 
are some of the issues discussed in this paper. 

We start with a quick resume of the boundary action method. We then generalize the 
method to arbitrary boundary conditions in section 3. In section 4, we apply this to the case 
of two parallel plates with Robin boundary conditions, one of them being a semi-infinite plate 
so that edge effects and diffraction at the edge become important. The analytical part of the 
calculation is explained in some detail, followed by the numerical evaluation of some of the 
integrals needed. A variety of values of the Robin parameter are considered. The case of 
two full plates (i.e., with no edge effects) with Robin boundary conditions [HI [12] and the 



geometry of a plate and half-plate but restricted to Neumann-Neumann conditions [HI [10] have 
been analyzed by other methods previously. We give a comparison of our results with these in 
the discussion section. Some other interesting features of our calculation are also pointed out 
there. The paper concludes with an appendix where some of the mathematical calculations we 
have done are outlined. 



2 Resume of boundary action method 

We start with a brief summary of the approach developed in O [H [5]. We will consider a 
scalar field theory with action S{(j)) in a cubic box of volume V, with y — )• oo eventually. 
For simplicity, we start off with a free field theory with a kinetic energy operator given by 
the Laplacian, i.e., S{(j)) = ^ f{d(j))'^. The quantity of interest for the Casimir energy is the 
functional integral (or partition function) 

Z = j[d4>] exp[-5(</))] (1) 

Our basic strategy is to consider the box as divided into a left region Vl and a right region 
Vr as shown in Fig.lll with an interface (shown as the dashed line) at, say, xi = b. Then 
the functional integral can be done in two stages. We integrate over all field configurations 
in Vl with a fixed value of the field (say, = 93) on the interface and similarly for the right 
region Vr. Each integration leads to a functional of (^, which we will denote by ^ l and ^ r 
respectively, up to constant factors Zi, Zr, so that 

Z = ZlZr [d(f] exp [-5^(99)] 
^L^R = exp[-SBiv)] (2) 

The final integration over the values of the field (p at the interface completes the functional 
integration in Q. If there is a plate placed at xi = b with the field vanishing on it, then 
this last integration is trivial; we just set ip = 0. If we have a plate with an opening, then 
ip = everywhere on the interface except at the opening. The integration of exp(— 55(93) ) 
with c/3 restricted to being nonzero only at the opening gives a boundary contribution to the 
partition function, which will include all the diffractive effects due to the opening. This is 
our approach in a nutshell. It is clear that the formalism can be extended to boundaries with 
many components with different configurations of apertures, edges, etc. There will be many 
contributory terms to Sb{^) and these would capture the various boundary effects. 

This stage-wise functional integration is easily implemented by writing (pL = (f'OL + Vl, 
where 0ol is a specific field configuration in Vl with the boundary value (poL -^ (f as xi ^ b. 
Explicitly, we may take it to be 

'AoL(a;) = / if{x') n-d^iGL{x,x') (3) 

JdVL 

where Gl{x,x') is the Green's function for D with Dirichlet boundary conditions, i.e., 
OGl{x,x') = d^^'{x — x'), Gl = if X or x' S OVl- In this equation, n ■ d denotes the 
derivative normal to the boundary OVl- The functional freedom for the fields in Vl is in r]L 



which is taken to vanish at all boundaries (including the apertures), so that the full field ^l 
has the required boundary behavior, (pi — )■ 99. This is consistent with the Dirichlet boundary 
condition for the Green's function Gl, which is also the propagator for rji. The functional 
integration over rji gives the partition function for the fields in Vl^ namely, Z^, evaluated with 
the Dirichlet boundary condition, and also ^l, which carries the (/^-dependence. A similar 
splitting 0R = (pQR + VR can be used for the right region Vr. The boundary action is then 
found to be 



Sb{v) 



ip{x) n- dx n- d^i [Gl{x, x') + Gr{x, x')] ip{x) 



(4) 



where the integration is over the common interface and where Gr also satisfies Dirichlet bound- 
ary conditions. 

We see that the value of the field at the interface, namely, (p is what is left free until the 
last integration. It is easy to impose a vanishing condition on this field wherever there is 
material plate at the interface. However, imposing a Neumann condition (or anything more 
general) is not straightforward. We may use a mode expansion of the fields in Vl in terms of 
eigenmodes of □, obeying more general boundary conditions, but then extracting a boundary 
action (which can easily accommodate apertures) becomes awkward. For example, if we want 
to use Neumann condition, we have to use bulk modes which obey this condition for the plate 
parts of the interface, but must have an arbitrary value ip on the openings, left free until the 
final integration. Finding such modes is a rather involved mathematical problem. We would 
like a formalism which has the facility of dealing with apertures and diffractive effects as easily 
as we can do for the Dirichlet case, yet accommodate general boundary conditions directly 
in terms of the boundary action, so that the bulk modes are easy to construct. In the next 
section, we will address this question. 

Two remarks are useful at this point. Although we used a free scalar field theory to 
illustrate the set-up, it is clear that we can use a similar formalism for an interacting field 
theory as well. Secondly, if we think of the xi-direction as time (with a Minkowski signature), 
the result of doing the functional integral with the specified value of the field, namely, 93 at 
xi = b, would be a wave functional (for some state) of the field theory. We have emphasized 
this by using the notation ^^(99). (If we take the left end of the box to be at xi = — cx), then 
we get the ground state wave functional.) So what we have is really a Euclidean version of 



Vl 


Vr 



Xl 



Figure 1: Schematics of splitting the functional integration into that over fields in two regions 
Vl and Vr 



this wave functional set-up. As shown in [3j, this wave functional can be obtained directly 
from the quantum effective action r[x]. We solve {5Tl/5x) = subject to the boundary value 
X ^ 93 as xi — )• 6. Then ^l = exp(— F^) where F^ is given by F[x] evaluated on this solution. 
(Obviously, similar statements hold for ^/j as well.) This gives another way to think of our 
procedure in a fully interacting quantum field theory. 

3 General boundary conditions 

We start again with the simple case of a free scalar field and consider the von Neumann theory 
of self-adjoint extensions [13]. This has recently been rephrased nicely by Asorey, Ibort and 
Marmo, so we can use their approach [TJj. For the Laplace operator in some region V, we 
consider the field (p and its normal derivative dn^- On the boundary, we consider square 
integrable functions, i.e., they form a Hilbert space. Thus the combination ip + idnP (where 
dn = n ■ d is the normal derivative) may be viewed as an element of this Hilbert space. The 
most general boundary condition, according to the von Neumann theory, is 

if + idn(p = U {ip - idnp) 



{if + idnP>){x) = (f)U{x,y){ip-idn(p){y) (5) 

where U is a unitary operator on the boundary Hilbert space, namely, on the space of L^- 
functions on the boundary, of which ip is an element. In the second line of the equation above, 
we emphasize this by writing U{x,y). Using the operator notation of the first line of (pi), we 
find 

dnP = -i i TTTT ] P = -K.p (6) 



u + i. 

/C is a hermitian operator and the transformation from {7 to /C is the so-called Cayley transform. 
There are two interesting limits, /C — t- and /C — t- 00. The first one, equivalent to U = 1, gives 
the Neumann boundary condition and the second one (U ^ —1) corresponds to the Dirichlet 
condition. (For the latter, it is better to divide ([g]) by /C and then take the limit.) These are 
special points. The case of /C being a constant (proportional to the identity on the Hilbert 
space) is the Robin condition. The general theorem is that, in the space of /C's, infinitesimally 
close to the two limits of Neumann and Dirichlet, the Laplace operator can have negative 
eigenvalues [Hj. 

The obvious question is whether this can affect our evaluation of the functional integral. 
To ensure consistency of the variational problem with the boundary conditions, extra surface 
terms may have to be added to the action, and one can ask if they are sufficient to avoid any 
pathologies. The answer, in general, is that the negative eigenvalues can lead to pathologies; 
but for certain types of /C, or range of eigenvalues for the same, we can have a stable situation. 
This is further commented upon in the discussion section. 

There is a different way to look at this problem, which also suggests the solution. Let us 
go back again to (pi) writing the partition function as 

Z = [det(-nL)]-2 [det(-nfl)]-5 [ [dp] ^ l{p) ^ r{p) 



= ZlZrZb (7) 



If the partition between Vl and Vr is real, like a plate, may be with openings, we express 
that as the vanishing of the field (j) = ip in the plate region of the partition, but free on the 
open region. Equivalently, we may say that the plate is represented by S{ip) where the delta 
function is only on the plate region. Explicitly, we can expand cp in term of modes on the open 
regions of the interface, and on the plate parts. The coefficients of the latter vanish by the 
delta function. The result then agrees with what we did in |3]. We may interpret this delta 
function as the operator representing the insertion of the plate, the plate operator. So the final 
integral, apart from the determinants, looks like 



Zb 



[dip] ^^li^) Spiatei^) ^R{ip) 



(8) 



Now consider what happens when we have a normal derivative. The normal derivative acts as 
a functional derivative on the ^'s. This can be seen easily from slicing up the path integral 
along the direction normal to the interface, the xi direction. The Euclidean action then looks 
like 



S = SiWi}) 



1 



d^x^ 



{ipN - ifN-lY , {^N-1 - ^N-2Y 



XN - XN-1 



XN-1 - XN-2 



+ (v^^)^ 



ipN is the boundary value ip. The boundary action is given by 

JV-l 



'Sb{v) 



/ Jl dipi exp(-S'({c^j}) 



(9) 



(10) 



Differentiating e with respect to (fN, we find 



6ip 



N 



XN - XN-l 



(11) 



So, going back to ([7||8]), we see that imposing a boundary condition ip' + Idp = on the plate 
is equivalent to imposing 

r r 

-— +ICip ^{ip) = (12) 

dip J 

To see how this requirement can be obtained, consider the integral 

I 
We can then write 

[dip] exp ( -- / iplCip] (dnip) "^ 



([dip] expl--! iplCipj ^ip] 



(13) 



/ [dip] exp 
[dp] exp 



ipKip 



5_ 

5ip 



m 



ipiCip] i-icp)np] (14) 



1 

2 
1 
2 

where, in the second line, we have done a partial integration. This equation can be rewritten 
as 

[d^]exp{-^ fiplCip] {dn^ + )Cip)^[^] = (15) 



Thus we do obtain the required vanishing of dnf + Jdp, showing that the operator representing 
the plate must be taken as exp(— ^ J cplCip). The boundary condition apphes only for the fields 
on the plate part of the boundary, so we must restrict the fields in J ipldp to be only on the 
plate part. Further, we can choose the boundary conditions independently on the left and 
right sides of the same plate, which have independent boundary fields tpi and ipR, so that 
the insertion of exp(— ^ f (p)C(p) is to be done separately for each. (This issue did not arise in 
previous calculations because we used Dirichlet conditions setting 99 = for both sides of the 
plate.) The fields on the aperture part are the same on both sides of the boundary. 

So the calculational algorithm is: Use Dirichlet conditions with left and right regions and 
fields (pL and (/)/j. These go to the boundary values ipL and ipR. Then calculate 

J \ ^ Jplate,L ^ Jplate,R / 

X-^[ifL]-^[fR] (16) 

We integrate over all fields (pi and ipR, but the (5- function ensures that the fields are the same 
on both sides of the aperture and the integrals in the exponents are restricted to the fields on 
the plate part of the boundary on the left and right sides. This formulation, at least partially, 
avoids the problem of potential negative eigenvalues since the bulk determinants leading to Zl 
and Zfi are always calculated with Dirichlet conditions. 

For the free theory, ^{(f) is of the form 

^((/7) =exp(-i f ^{x)M{x,x')^{x')] (17) 

where M{x,x') = n- dxn ■ dx' G{x,x'). Thus, in this case, the expression ( 16 ) for Zb becomes 

Zb = [d^L dipn] Saperturei^L " ^r) e~^^ 

Sb = 7, V^lICl^l + i; ^rK,r^r + -1 (pMifi (18) 



^ J plate, L -^ J plate, R ^ J boundary 

Once again, the integrals involving /C's are only over the plate regions of the boundary, while 
the last integral is over all boundary. If the boundary has disconnected components, as would 



be the case for, say, parallel plates, then all such components must be included in (18). We 
are now in a position to apply this method of calculation to a specific example. 

4 Plate and Half-plate 

4.1 Modes on the boundary and the boundary action 

We will consider the arrangement shown in Fig. [2] where we have an infinite plate and a semi- 
infinite plate parallel to it, separated in the 2;i-direction by a distance h. The case when the 
field has Dirichlet boundary conditions with </> = on the plates was analyzed previously [3]. 
Here we will consider more general boundary conditions. The bulk contribution from the left 
and right regions, namely, Zl Zr, which in the present case is also calculated with Dirichlet 
boundary conditions, is the same as before. For the boundary contributions, first of all, we 



need the boundary action. The Dirichlet propagator in the left region, between the two plates, 
is given by 



G{x, x) 



(2^ 



JH-^-'^)Gd{uj,xi,x\) 



- N{u) \e{xi - x[) (e^(^i-<) + e^^ e-^^^^-^'i)) 



ix)(x\^x'^) _|_ 26w -w{xx^x'^ 



where Q{x\ — x'^ is the step function and 



o2bLC _ ]^ 



(19) 



(20) 



Further, uj = \fW = y^/cg + /c^ + fe 



3' 



r,T 



{xq,X2.,x^). We will take the open part of the 



second plate to be in the range X2 > 0. Calculating the normal derivatives, we find 



1 



Lpi{x) {uj coth.hu} + K,i)^y Lpj{x')+ / ipi{x){-uj csch.hu}) ^^^, (fn{x') 
+ 7: I fii{x) {ujcothhu} + ICL)^^f (fjj{x') + - / fR{x) {u + K,r)^,^, lpr{x) 



Sb 



(21) 

where (pi refers to the field on plate I, at xi = 0, and ipn to that on the left side {Vl side) of 
the interface at xi = 6 (which includes the half-plate and aperture) and ipR to the field again 
at the same interface, but on the right side (the Vr side). 

We will now need to specify the operators /C/, K-l and K.r. For illustrating the techniques 
outlined, we will choose a simple Robin boundary condition /C/ = nj = constant for the left 
plate (plate I). The special cases k/ = and k/ — )■ oo will correspond to the Neumann and 
Dirichlet limits for the left plate. To specify ICl and ICr, it is useful to separate ipji into modes 




Figure 2: The arrangement of a Robin plate and a Robin half-plate with values of fields 
indicated 



which have support on the aperture, modes with support on the plate, and the value of the 
field at the edge of plate II, X2 = 0. This is done by writing 



y^ii 



(2^ 



ik-x 



dp 



2p 



IT p^ 

2 



+ 



fc2 + fc2 



{k2- 



le 



■ak,p) 



dp 



2p 



TT p"^ — {k2 + ie)2 



<p2Lik,p) 



(22) 



Here k = {ko, k^), the two directions transverse to the edge of the plate, as well as transverse 
to the plate itself. We define the fields 



')2L{X) 

e(^) 

p{x) 



(2^ 
d^k 



2^e 

TT 

vr 



ik-x 



ik-x 



sm{px2) <l)2Lik,p) 



sin(px2) i{k,p) 



(2vr)- JO 

(2vr)2" ^^^M el^>^ X2<0 



(23) 



The field ^ is the field at the aperture (x2 > 0) but vanishing at the edge X2 = 0, (j)2L is the 
field on the plate {x2 < 0), vanishing at 2:2 = and p is essentially the value of the field at 
X2 = 0, but continued in a very specific way, with no additional functional degrees of freedom, 
to X2 > and X2 < 0. (While the value of the field at X2 = ±00 will not be relevant for us, 
the value of the field vanishing at X2 = ±c« will make some of the integrals easier and better 
defined. The continuation of the field value at the edge, namely p, to other values of X2 has 
been done in one particular way which ensures this.) 



Carrying out the p-integration in ( 22 ) shows that 

^{x) + p{x) 
(pn{x) = { 4>2l{x) + p{x) 

p(xo,X3,X2 = 0) 



X2 > 
X2 < 

3:2 = 



(24) 



The amplitudes of the modes in (22), namely, C{k,p), p{k), and 4'2Lik,p) constitute the func- 



tional freedom in the value of the field at the boundary. In a similar way, we can write 

Vr{x) = { (p2R{x) + p{x) X2<0 (25) 



({x) + p{x) 

4'2r{x) + p{x) 
p{xo,X3,X2 = 0) 



X2 > 
X2 < 
X2 = 



where 



52R 



d^k 



2^P ik-x 

vr 



sin(px2) 02i?(fc) 



For if I, since there is a single plate for all X2, we can use a simple mode expansion, 

d^k 



^i 



ik-x 



(27r)3 



lik) 



(26) 



(27) 



The normalization for sin(px2) is 

/ dx2 sin(px2) sm{qx2) = -6{p - q) (28) 

Jo ^ 



One advantage in parametrizing the fields as in (22) is that this already takes care of the 
5- function 6aperture{^ii — ^r) enforcing equality of fields on the aperture X2 > 0. 

We may now regard ICl as an operator on the fields {(p2L, p) and likewise K-r as an operator 
on {4)2Ri p)- While ^ and 4)21 (and ^ and 4)2r) are orthogonal to each other, p is not orthogonal 
to ^ or (f)2L (or 4>2r)- Nevertheless, (^, (j)2L, p) and (^, (/)2_r, p) form a complete basis for fields at 
xi = b. There is clearly an infinity of choices possible for ICl and ICr, but, once again, for our 
analysis, we will make a simple choice. We will take ICl = kl = constant for all the modes (J)2l 
and /C/j = kr = constant for all the modes (/>2ij. In other words, the operators ICl and ICji are 
diagonal with the same eigenvalue kl (respectively kr) for all (j!)2l (respectively (J)2r)- This is 
almost like a Robin condition for the plate on the right as well. We say "almost" because the 
situation with p is a bit tricky. It corresponds to the field just at the edge X2 = 0. In principle, 
the value of ICl for p can be different from the values of the same operator for (/>2l , even after 
we have chosen the latter to be the same for all modes <^2L- A similar statement applies to 
/Cij for p versus ICr for (j)2R- In addition, we can have nondiagonal terms mixing p with (J)2l 
and (p2R- Thus even with the simplifying choices we have made, there are many parameters 



specifying the boundary behaviour. Summarizing, the boundary action (18) for this geometry, 
with the choices we have made, is 

Sb = - (pi{x){u}cothbuj + Ki)^^^'(fj{x') + / (pi{x){-ujcschbu})^y(pii{x') 
+ 2 / V//(a;) {u)cothbuj)^yipij{x') + - ^r{x) {uj)x,x'Vr{x') 

+ 2'^L / (f>2L4>2L + ^Kr / (f>2R4>2R+ 2*^ j P P 

+C2L / (t>2LP + C2R / (t>2R P (29) 



4.2 Simplification of the boundary action 

We are now in a position to evaluate the action in terms of the mode expansions we have given 
and then carry out the functional integrals. We will only give the final form of the boundary 
action here, relegating the details to the appendix. 

We will separate the terms into two categories. Because p{k) depends only on two dimen- 
sions, it is easier to integrate out the other fields first and leave the p integral to the end. 
Towards this, we shall first simplify the terms involving ipj, (j)2L, 4'2R. and ^. In this way, we 



separate out the terms in (29) as 



Sb = S^i'> + 5g^ 



- / (/'2L A^2L <A2L + 2 / 4'2rM2R(I)2R+ / 4'2lQ2L^C+ / '/'2R A^2R^ ? (30) 



+ 2 



10 



S'i^ = - pMpP+ h)lNipP+ <I)2lQ2LpP+ (t>2RQ2RpP+ pQp^i (31) 

The various coefficient functions, sucli as Mi, Ni2l^ etc., arise naturally from the restriction 
of the coefficient functions in (29) to the appropriate modes. We do not give their expressions 
here; they are given in the appendix. Integrating over </>i, (J)2l and (J)2r yields 

- log Zb = -Tr log(L^ coth hoo + ki) + -Tr log(Mi) + -Tr \og{MR) - log Z'b 



(32) 



Z'b 
S'b 



[d^dp] e 



-s' 



1 



lk,p,q 



a-k,p)M^{p,q)C{k,q) 



+ I p{-k)Np^{k,p)i{k,p) + -j^ p{-k)Mpp{k) 



k,p 



' k,p,q 



k,p,q 



N2Ldk, P, P) a-k, P) - p{-k) N2Lp{k, p) 



N2Ldk, q, q) C{k, q) - p{k) N2Lp{K q) 



N2R^{k, p, p) i{-k, p) + p{-k) N2Rp{k, p) 



NiR^k, q, q') i{k, q') + p{k) N2Rp{k, q) 



G2Lik,p,> 



G2R{k,p,, 



The indicated integrations are done with the measures 

r _ r d^k dpdq r _ 

Jk,p,q J (2vr)2 TT TT ' J p> 



'(V 
^ 



etc. 



The other quantities in (32, 33) are 

Ml = 2TT{Hp + KL)6{p-q) + {AmL)pq 



(33) 



(34) 



(35) 



(AMl) 



pg 



4pq 



ds 



vr 



H. 



+ 



i?. 



+ 



H„ 



{p^ — s'^){q^ — s^) [s^ — p'^){q^ — p"^) {s^ — q'^){p'^ — q^) 



Mr 

(AMR)pg 



2^ {LOp + kr) 5{p -q) + (AMR)pg 
'^^ ds 



(36) 



Apq 







vr 



W.c 



+ 



UJr, 



+ 



iOn 



{p^ — s'^){q^ — s"^) {s^ — p'^){q'^ — p^) {s"^ — q^){p^ — q^) 



27r {Hp + ojp)6{p -q) + {AM^)pq 



(37) 



(AMg) 



pq 



4pq 



Mk\' 



ds 



^ 



{Hs + U}s 



+ 



[Hp + ujp) 



+ 



{Hq + i^q) 



[p^ — s'^){q^ — s^) (.s^ — p'^){q^ — p^) (s^ — q'^){p'^ — q^) 



ds Hs + ujs 

2^ w! 



+ 



Co 



(38) 
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In these formulae, 



Hr,=UJr,[ ^^ I , UJ = k +p 



(39) 



^ ^ ^^UpCothbujp + Kj ^ 

The quantities G2L and G2R are the Green's functions or inverses of the kernels M.l and M/j 
respectively. Explicitly, they are of the form 



G2l{p,< 



G2Rip,i 



TT 5{p - q) 
dp' 



1 



1 



-(AMl)p, 



1 



1 



2{Hp + kl) 2{Hp + kl) 
f°° dp' 1 

2(a;p + kr) ~ 2{L0p + kr) ^^^^^pi 

°° dn' 1 1 

+ / —^, : .(AMr),p, 



2{Hg + kl) 



(AMl) 



p'g 



2{Hg + kl: 



+ •••(40) 



vr 6{p - q) 
dp' 



2{ujq + Kr) 



Iq TT 2{u}p + KR 

Finally, we also have 

Npl:{k,p) = 4\k\p 

N2Ldk,P,P) = -(AMi)p, 

N2R^{k,p,p') = -{AMR)pg 



'=^ ds 
^ 



2{u}pf + kr) 

Hs + tJs) (i:/"p + Up) 



(AM 



i?jp'g 



2(a;q + kr] 



+ •••(41) 



w^ 



p2 _ g2 



4pp' 



(^p/2 _ p2'j 



{p'^ + k^)I{k,p') - {q^ + k^)I{k,, 



(42) 



where 



N2Lp{k,p) = Mk\p 
N2Rp{k,p) = -A\k\p 

I{k,p) = 



ds / Hs Hr, 



TT V^s ^pj P^ - S^ 



1 C2L 



2|A:|a;2 



X(/c,p) + 



C2R 



f^2 _i_ p2 _|_ p 



2'Kpy k'^ + p^ 



1^2 _i_ p2 _ p ^ 



(43) 



(44) 



The equalities N2l^ = — AM^ and N2R£_ = — AMr show that N2l^, A^2i?$ are of the same 
order in diffractive effect as the AM's, an issue which is relevant when we do the expansions 
explained below. 



4.3 Evaluation of the Casimir energy 

We are now in a position to write down the free energy, which will have several terms corre- 
sponding to the contributions from the left and right bulk regions and integrations over (pi, 
4>2L-, <t>2Ri etc. First, in extracting the free energy, we note that the /c-dependence is the same 
for all terms and that the overall integration factor representing the trace is LQL^{d^k/{2iT)'^). 
With — log Zb = Lq F, the free energy is given by 



F — Fbuik + Fj + F2L + F2R + F^ + Fp 



(45) 
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We have added in Fbuik, which is the contribution from the bulk determinants (from Z^, Zji), 



calculated with Dirichlet conditions. The remaining terms in (45) arise from the boundary 



action and correspond to the terms on the right hand side of ( 32 ) or the result of integrating over 



the terms on the right hand side of ( 33 ) . All these contributions have divergences corresponding 
to free space with no plates. They can be identified as the 6 — )• oo limit of the expressions here. 
Thus, in the following calculations, the renormalization of all terms will be done by subtracting 
the 5 — )> oo limit. 



A note on the expansion scheme 

The exact calculation of the relevant determinants and the free energy is still very involved 
and an expansion scheme is needed to get a good approximation for several of the terms in 
(45). In our earlier work [3], we noticed that the relevant operator, say M, had the structure 
M.^^'6{p — q) + AMpg where M^^^ is diagonal as indicated. It was referred to as a pole 



^pq 



term and AM was designated a cut-term, based on the method of calculation we had used 
[3]. More appropriately, M*^*^) was a "direct term" giving the geometrical optics approximation 
while the "diffractive term" AM captured the effects of diffraction [^ . We then carried out an 
expansion in powers of the diffractive term, i.e.. 



Trlog(M(°) + AM) = TrlogM^") + Tr((M(°))"^ AM) 

--Tr((M(°))"^ AM (M(°))~^ AM) + • 



(46) 



Although there is really no parameter controlling this expansion, numerically the higher and 



higher order diffractive contributions in (46) were smaller and smaller, and a sensible truncation 
was possible. 



We propose to do a similar expansion here. The propagators in (40, 41), as we have 
presented them, naturally show this expansion. Likewise, N2l^ = —AM/,, N2r^ = — AM/j are 
to be considered as of the first order in the diffractive effect. With this understanding, we can 



now proceed to the individual terms in (45). 



The bulk contribution 

The bulk contribution to the free energy is 

1 



nuik = -^Ls{w^ + w2)j ^2^^3 



d^k 



log 1 



-2bui 



(47) 



where Wi is the width of the aperture, W2 is the width of the half-plate on the right, with 

L2 = Wi+W2. 



The contribution from fields on the left plate (plate I) 

For Fj which is the contribution from the left plate, we have 

Fj = -Trlog((^fc coth6(^A; + nj) — (6 — ;■ 00 limit) 
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MWI + W2) 



(2^ 



log 



LOk coth buk + Hi 

Uk + Hi 



(48) 



We have also carried out the renormalization by subtracting the 6 — )• 00 limit. 



The contribution from fields on the left side of the half-plate 

The term F2L corresponds to oTrlogM^, and can be simplified as 



1 . ... /" d^k dp 

— (6 — )• 00 limit) 



1 f d^k 

log {Hp + kl) + -L3 J --pTr log 



1 + 



2iHp + kl) 



AMl 



(49) 



where we have expanded the second term in a series in AM^,. In expanding the last term, we 
get the trace of products of {l/2{Hp + kl)) AM/,. The "Tr" stands for setting the labels p = q 
for the two end terms and integrating with dp/n. The first few terms in the expansion are 



2L = 0-^3^2 / TTT^— log , 

2 J (27r)2 vr \u}p + kl 



.(i) _ L3 f d^k dp 



2L 



7(2) 

2L 



1 



2 J (27r)2 vr 2{Hp + kl\ 
L3 f d'^k dp dq 1 



(AMl)pp - (5 -^ 00 limit) 



(50) 

(51) 



4 J (27r)2 vr vr 2{Hp + kl) 



(AM 



Ljpg 



2{Hq + kl; 



(AML)gp - (6 ^00 limit) (52) 



The contribution from fields on the right side of the half-plate 

The term F2R corresponds to ^TrlogM/j; it has no 6-dependence and hence it is eliminated 
by renormalization. 



The contribution from fields on the aperture 



The next term of interest is the contribution from the ^-integration. The ^^ terms in (33) can 
be collected together as 



S> 



1 



ss 



k,p,q 



a-k,p) 



liP^Q) 



N2Li(.k,p ,p) G2L{k,p' ,q) N2L^{k,q' ,q) 



N2Ri^{k,p,p)G2R{k,p,q)N2B^{k,q,q) 



p',q' 



ak,q) 



a-k,p)2{Hp + u;p)ak,p)+ I a-k,p)iAM^)pgak,q) 

k,p Jk,p,q 



(53) 



where 



(AMg) 



pq 



(AM^) 



pq 



N2Li{k,p ,p) G2L{k,p ,q) N2L^{k,q ,q) 



p',q' 
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N2R^{k,p,p) G2R{k,p, q) N2R^{k, q, q) 



t(i)^ 



(AM4J = - [ N2Ldk,P,p)w7TA. .N2Ldk,P,q) 



(AM^)^ 






N2R^{k,p,q) 



(54) 
(55) 



(56) 



For the contribution to free energy upon integration over the ^'s we find 



I ^ ^^l- (Pk dp ^ / .. N 1 . [ (fk ^ ^ , 



(27r)2 IT 

— (6 — )• oo hniit) 
We can simphfy this further by noting that 



Hp + iOp 



(27r)2 



2{Hp + iOp 



-(AM^)) 



{Hp + ujp)i,^ 
The term F^ thus sphts up as 

1 



1 



1-e 



— 2buJr, 



UJp + K/ 



Wp coth hojp + K/ 



(57) 



(58) 



F^ = Fl"" + F,^^^ + F'/> + ■■■ 






F. 



.(0) . 

? --2--^V (^3 

(1) L3 I' d?k dp 1 



1 



-L3W1 



d^k 



5 2 7 (27r)2 vr 2(i7p + Wp 



2'-^ ^ J (27r)3 
(AM^)pp - (6 -^ 00 hmit) 



log 



bJk coth 6cjfc + Kl 
UJk + Kl 



(59) 
(60) 



F 



(2) -^3 /" c?^^ dp dq 

1 J (27r)2 W 

L3 /■ (i^A; dp 



1 



2(i?p + u)p 
1 



■(AM^) 



1 



-(AM) 



2 7 (2^)2 vr 2(i7p + ujp) ^^^« V 



2(//p + LOp) 

— (6 —7- 00 limit) 



qp 



(61) 



Notice that the terms in (59) cancel against similar terms in F[)uik and -F/. 



The contribution from the field at the edge (field p) 

The contributions we have evaluated so far correspond to the fields (j)2L^ 4'2R and ^ which 
vanish at the edge of the half-plate. If the parameter cq — t- 00, these are the only terms we 
have. For finite values of cq, we have the contribution from p as well. To evaluate this term, we 
need the kernel for the pp-teim in the boundary action. In addition to the pp-ieim manifestly 



displayed in ( 33 ) , there will be additional terms from the integration over ^ because of the ^-p 



terms in (33). For this, first define the propagator for ^ given by (53) 



G.{k) = n6ip,q) 



1 



1 



2{Hp + ujp) 



2[Hp + ujp 



■(AM^) 



1 



pq 



2iH„ + UJa 
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f°° dn' I ~ 1 ~ 1 

The integration over the ^'s thus leads to the p-dependent terms 

Sp = l Lpi-^) fMp(fc) + AM J pik) (63) 

where 

AMp = -/ N2Lp{k,p)G2L{k,p,q)N2Lp{k,q)- N2Rp{k,p)G2R{k,p,q)N2Rp{k,q) 

Jp,q Jp,q 



N,^i-k,p)G^{k,p,q)N,^{k,q) (64) 

'p,g 

Npi = Np^- f N2Ldk,p',P)G2L{k,p',q')N2Lp{k,q) 

Jp',q' 

N2Rdk,p' ,p)G2R{k,p' ,q')N2Rp{k,q') (65) 



p ,q 

The quantity AMp naturally has an expansion which follows from the expansion of the prop- 
agators G2L, G2R and G^. Notice that we start off with at least two powers of N2LS,-: -^2_R$ or 
A'^p^, so that there is no first order term. The expansion for AMp is thus 

AMp = AM{,2) + AmJ,^) + • • • 
AM(2) = - I N2Lp{k,p) \ . N2Lp{k,p) - I N2Rp{k,p)—^ ,N2Rp{k,p) 

'^ Jp 2{Hp + KL) Jp ^(Wp + Kr) 

^pd-k,P)w7TF^—-xNpdk,p) (66) 

Z[np + LOp) 

All terms in these expressions depend on k; we have not indicated this explicitly in the argu- 
ments of various functions to avoid too much clutter. Once the integrations over Pj_g^ etc. are 



done, Mp + AMp is only a function of k, so that the free energy resulting from (63) can be 

written as 

Li r cPk 

Yj (2^ 
Thus approximations to this term will be of the form 



^P = ^ / A^^ l°g (Mp + AMp) (67) 



i^r = ^ I ^ log (Mp) - (6 ^ 00 limit) 



FW+i^W = ^y"^log(Mp + AM(2))-(6^oolimit) (68) 



4.4 Numerical estimates 



The terms proportional to Wi, the width of the aperture, cancel out between (47), (48) and 



(59). Thus we are left with bulk terms which only involve the facing area L3W2 between the 
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two plates. This is given by part of Fbuik in (47), part of Fj in (48) and ^2^^ in (|50|) and is 
equal to 



LsVFo r d^k 



(2^) 



log ( 1 - e-^'"'^ 



+ log \ + log 



UJk + Kl 



^k + I^L 



(69) 



There are also many edge-dependent terms to consider. The fields ^, (J)2l vanish at the edge 
X2 = 0. For an aperture and a half-plate of finite widths Wi and W2-, we should use discrete 
modes to evaluate the free energy, with a sum replacing the integral over p or k2- In the limit 
of large Wi, W2, the sum over the discrete modes can be approximated by an integral, as in an 



Euler-Maclaurin summation formula, which gives the area-dependent term in (69). But there 



is also a subdominant edge-dependent term in the summation formula [3j; we will denote this 
term by Ff^]^- 



The remaining edge-dependent terms in the energy are given by Fg^, + F2L ~^ 



from 



(|5lD, (|52|),etc., Fi^) + Ff^ + 



p lie I from (60), (61), etc., and Fp from (68). All these quantities 

depen3~on the parameters k/, kl-, hr, cq, C2l and C2r. We will now proceed to the numerical 
estimate of these integrals for some choices of these parameters. 



The area-dependent contribution 



The area-dependent part of the Casimir energy (69), from scaling out b, is of the form 

L3W2 vr2 



63 1440 



£area{bKl ,bKL) 



(70) 



We will consider Dirichlet condition for the left plate (k/ — >• cxd), Robin condition for the 
right (with finite bn^ = m), or the other way {kl — )• 00, bnj = m finite), and then Robin 
condition for both (with bKj = bui = fn). The values for these cases are given in Table IT} The 
Dirichlet-Robin and Robin-Dirichlet cases give the same area-dependent term for the energy. 
The Robin-Robin case with tti = corresponds to Neumann condition on both plates. The 



Table 1: The area-dependent contributions Sarea for different boundary conditions 



m 





1 


2 


3 


4 


5 


£rd{oc, m) = SoRim, 00) 
£RR{m,m) 


0.875 
-1 


0.092 
-0.089 


-0.212 
0.115 


-0.381 
-0.200 


-0.489 
-0.283 


-0.564 
-0.354 


m 


6 


7 


8 


9 


10 


00 


£rd{oo, m) = SoRim, 00) 
£RR{m,m) 


-0.620 
-0.414 


-0.663 

-0.464 


-0.698 
-0.507 


-0.725 
-0.544 


-0.749 
-0.575 


-1 
-1 
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£dd(oo,oo) 



Figure 3: The area-dependent contribution to energy, Sarea, for different boundary conditions, 
as a function of m, showing asymptotic approach to the Dirichlet-Dirichlet case. 



limit of m — )■ oo should correspond to the case with Dirichlet conditions on both plates which 
is, of course, well known. As expected, and as Fig. [3] shows, the values found here approach 
that value {£ = —1) asymptotically. 



The edge-contributions (^vithout p) 



(0) 



The edge-dependent contributions from F2L and F^ (including F^j^ ) are of the form 



L 



Fedge = p" £edge{bKl ,bKL,bKR) 



(71) 



> ooj 
bKR 



Dirichlet- 
finite 



m= 



We consider the cases of Robin-Dirichlet (6k/ = m = finite value, k/,,kr 
Robin (k/ — ?■ 00, Bkl = bKFt= m = finite value) and Robin-Robin (6k/ = bKL 
value) . The values Sedge for a range of m are given in Tables [2I [sl [4| 

We have calculated these to the second order in the expansion of the Tr log terms, as in 
(50 52), (59 61). These do not include the contribution from the p- integration which will be 



given separately. We also graphically display the leading order term, the first order term, the 
second order term and the total edge contribution to this order for the three cases in Figs.|4| 
[5| [6| The dashed line is the total contribution, up to second order, and the dotted line is the 
Dirichlet-Dirichlet limit. 
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Table 2: The edge contribution to the Casimir energy, Robin-Dirichlet case, without Fp 



m 





1 


2 


3 


4 


5 


^eV,fio("^'°°'°°) 


0.00448 


-0.00079 


-0.00241 


-0.00325 


-0.00377 


-0.00412 


^eW,fio("^'°°'°°) 


-0.00030 


0.00005 


0.00016 


0.00022 


0.00025 


0.00028 


^ed3e,fio("^'°°'°°) 


0.00180 


0.00121 


0.00100 


0.00088 


0.00081 


0.00075 


^idgeM^^^^^) 


0.00249 


-0.00033 


-0.00112 


-0.00152 


-0.00176 


-0.00193 


m 


6 


7 


8 


9 


10 


oo 


^edge,RDi°^^^^^) 


-0.00439 


-0.00456 


-0.00471 


-0.00483 


-0.00493 


-0.00598 


^edge,RDi^^^^^) 


0.00203 


0.00212 


0.00219 


0.00224 


0.00229 


0.00277 


£edge,RDi^^^^^) 


0.00029 


0.00031 


0.00032 


0.00032 


0.00033 


0.00040 


SidgtRDi^^^^^) 


-0.00205 


-0.00214 


-0.00221 


-0.00227 


-0.00232 


-0.00280 



0.004 



0.002 



-0.002 - 



-0.004 




Figure 4: The edge-dependent contributions from F2L and F^ for Robin-Dirichlet case, 

Sedge,B,Di'm, 00, 00). 
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Table 3: The edge contribution to the Casimir energy, Dirichlet-Robin case, without Fp 



m 





1 


2 


3 


4 


5 


^edge,DRi^^^^^) 


-0.01644 


-0.01117 


-0.00954 


-0.00871 


-0.00819 


-0.00784 


^edge,DRi^^^^^) 


0.00790 


0.00594 


0.00516 


0.00471 


0.00441 


0.00419 


^edge,OR(°°'"^'"^') 


0.00180 


0.00121 


0.00100 


0.00088 


0.00081 


0.00075 


£XlDRioc,m,m) 


-0.00675 


-0.00402 


-0.00338 


-0.00312 


-0.00298 


-0.00290 


m 


6 


7 


8 


9 


10 


oo 


^edge,DRi°^^^^^) 


-0.00759 


-0.00740 


-0.00725 


-0.00712 


-0.00702 


-0.00598 


^edge,DRi^^^^^) 


0.00403 


0.00390 


0.00379 


0.00371 


0.00364 


0.00277 


^edgcDR^^^^^^) 


0.00071 


0.00068 


0.00065 


0.00063 


0.00061 


0.00040 


^iteM^^^M 


-0.00285 


-0.00282 


-0.00280 


-0.00278 


-0.00277 


-0.00280 



0.005 


^~~--~~,_____l-st order 
i.____^^^ 2-7id order 
























III 1 1 1 1 1 1 1 1 1 1 1 1 1 




2 


4 


6 




8 


10 


0.005 


-,*•*' Total 

^-'■^'^ 0-th order 


■ 








0.010 
0.015 











Figure 5: The edge-dependent contributions from F2L and F^ for Dirichlet-Robin case, 

Sedge,DRioO,m,m). 
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Table 4: The edge contribution to the Casimir energy, Robin-Robin case, without Fp 



m 





1 


2 


3 


4 


5 


^edge,RR(.^^ ^^ ^) 


0.01495 


-0.00084 


-0.00340 


-0.00441 


-0.00492 


-0.00521 


^edge,RRi^^ ^^ ^) 


-0.00515 


0.00058 


0.00187 


0.00238 


0.00263 


0.00277 


^!dge,RRi^' ^^ ^) 


-0.00159 


0.00006 


0.00034 


0.00043 


0.00047 


0.00049 


£iteM^^^^^) 


0.00820 


-0.00020 


-0.00120 


-0.00160 


-0.00182 


-0.00196 


m 


6 


7 


8 


9 


10 


oo 


£edge,RRi^^ "^' ^) 


-0.00540 


-0.00552 


-0.00561 


-0.00568 


-0.00573 


-0.00598 


£edge,RRi^^ ^^ ^) 


0.00285 


0.00289 


0.00292 


0.00294 


0.00295 


0.00277 


£ldge,RR{m, m, m) 


0.00049 


0.00050 


0.00050 


0.00050 


0.00050 


0.00040 


^iteM^^^^^) 


-0.00206 


-0.00213 


-0.00220 


-0.00224 


-0.00228 


-0.00280 



0.015 


















0.010 


















0.005 


1 \ 

1 \ 
\\ 
- \\ 






1-st order 
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'' 


2-nd order 












. 


1 
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x" 
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8 


10 


-0.005 


- 




-V, 


0-th ordeT' 


















Figure 6: The edge-dependent contributions from F2L and F^ for Robin-Robin case, 

£edge,RR{m,m,m). 
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The edge-contributions from p 

The edge contribution from p is zero for the case of cq — >• oo. Thus for the case of Robin 
conditions (or anything else) on the left plate and Dirichlet conditions on the half-plate, the 
results given so far suffice. We will consider the inclusion of the contribution from p for the 
Dirichlet-Robin case (k/ — )• cxd). For simplicity, we shall consider the case C2l = C2i? = 0. We 
will need to specify cq as well. We take it to be cq = {kl + i^r), since this is what would 
be naturally considered if we take p to be part of 02L and 02K- (There is no compelling 
reason for this choice; it is one case worth considering and easy enough to calculate.) Thus 
the contribution from p is of the form 



hi 
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£p{bKi,bKL,bK.R,bco) 



(72) 



We will calculate this to the second order as well, which means that we will evaluate the 
expression in (68). These values, for a range of parameters, are given in Table pi We also 
show the contribution from p as compared to the other edge contributions as a function of the 
parameter m in graphs Figs.lTl ISl 



Table 5: The edge contribution from p to the Casimir energy 



m 





1 


2 


3 


4 


5 


£p,DR{oo,m,m,2m) 
£p,RR{m,m,m,2m) 


0.00336 
-0.00553 


0.00356 
0.00185 


0.00311 
0.00102 


0.00271 
0.00127 


0.00239 
0.00133 


0.00214 
0.00132 


m 


6 


7 


8 


9 


10 


oo 


Sp,DR{oo,m,m,2m) 
£p,rrW^^ m,m, 2m) 


0.00193 
0.00128 


0.00177 
0.00123 


0.00162 
0.00118 


0.00150 
0.00113 


0.00140 
0.00108 







5 Discussion 

We have generalized our previous formulation for the calculation of diffractive effects in Casimir 
energy in the case of boundaries with edges and apertures to include general boundary con- 
ditions. As a specific example we have analyzed the geometry of two parallel plates, one of 
which is a semi-infinite so that there are edge effects and diffraction. We have considered a 
variety of boundary conditions and almost all the results are new, not calculated previously by 
any method. There are a few special cases for which results using other methods are known 
and a comparison with the literature is possible. 

The first case is where we have two full parallel plates. The Casimir energy for a massless 
scalar field subject to Robin boundary conditions on two infinitely long parallel plates has been 
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0.002 



-0.002 



-0.004 



-0.006 



From p 




Total without p 



Figure 7: The edge-dependent contributions, including p, for Dirichlet-Robin case, 
£edge,DR{oo, m, m, 2m). 



0.008 



0.006 



0.004 



0.002 



-0.002 



-0.004 




Frqrnj}_ 




Total without p 



Figure 8: The edge-dependent contributions, 
£edge,RRim,m,m,2m). 



including p for Robin-Robin case. 



analyzed before in [TTl [12] . Of particular interest is [12] , where a boundary term was added to 
impose Robin conditions at the level of Green's functions. The idea is similar to ours although 
our argument for such a term, which is based on the wave functional interpretation is different 
and somewhat more general. If we ignore the diffractive contributions in our calculation and 
replace the facing area of the plates by the full area our results are in agreement with |11[ IT2] . 
(The graph in [12] shows £/m^ rather than £ as in our graph.) 

The MIT group has analyzed, using methods of scattering theory, the edge contribution 
to the Casimir energy in a number of different geometries of boundaries with edges |9i llOj . 
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One particular case is exactly the geometry we have but restricted to DD and NN conditions. 
Our results for the DD, ND, DN and NN cases are summarized in Table [6j These values for 
the DD and NN cases are in good agreement with the results in [TJ [9l [10]. In particular, the 
numerical values obtained from the expressions given in |10] are ^£)^^ = —2.63 x 10^'^ and 
£%f^ = 2.97 X 10^^. In our previous paper [3] we calculated the edge contribution up to the 
5th order in diffraction, £'£,f^ = —2.68 x 10^^, and the agreement with [10] is even better than 
is apparent from Table [6j In making this comparison, it is important to keep in mind that the 
Neumann boundary condition corresponds to integrating over all boundary fields, including 
p (with Ki = hl = Kfi = cq = Qi)] hence the contribution from p must also be taken into 
account in comparing the values. We also note that the DD case was initially treated using 
numerical worldline methods in [7]. The value obtained was S'j^f) = —2.62 x 10~^; again the 
later calculations [3l|9l|10] are consistent with this result. 

Another interesting result which emerges from our calculations is that the diffractive effects 
are always of the opposite sign to the nondiffractive (and leading) contribution. This had been 
noticed already in our previous work O [5l Ej where we considered Dirichlet-Dirichlet boundary 
conditions. This continues to hold with the more general boundary conditions discussed here, 
as is strikingly clear from the graphs for the edge contributions. There is presumably a general 
reason for this, it is a point worth exploring. 

Finally, a short remark on the issue of negative eigenvalues: The bulk determinants in our 
approach are calculated with Dirichlet boundary conditions. There are no negative eigenvalues 
to worry about for this calculation. However, this issue is not totally eliminated; it has a 
lingering effect on the boundary action. Notice that the signs of k/, kl, etc. are important. 
These parameters occur in the boundary action and negative values for these can lead to 



instabilities. For example, the propagators in (40, 41) clearly display this possibility. We have 



to conclude that acceptable boundary conditions which encode the boundary effects of real 



Table 6: Summary for DD, ND, DN and NN cases 





DD 


ND 


DN 


NN 


'^area 

po 

edge 

Pl 

edge 

P2 
edqe 


-1 

-|g = -0.00598 
0.00277 
0.00040 


7 

IS? = 0.00448 
-0.00170 
-0.00030 


7 

-^iJ = -0.01645 
0.00790 
0.00180 


-1 

f|g = 0.01495 
-0.00515 
-0.00060 


ctotal 
edge 

^edge,p 


-0.00280 



0.0048 



-0.00675 
0.00336 


0.00820 
-0.00553 


^edge 

(Total) 


-0.00280 


0.00249 


-0.00339 


0.00267 
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material plates must be such that the operators K, have positive eigenvalues. One may still ask 
the question whether there is any meaningful physical interpretation for the mathematically 
acceptable case of negative eigenvalues. This will be discussed in a separate paper J15 j. 

This research was supported in part by NSF grants PHY-1068172 and PHY-1213380 and 
by PSC-CUNY grants. 

Appendix 

In this appendix, we work out some of the simplification of the terms in the boundary action. 



Using the mode expansions in (24, 25, 26) and (27), we find 
(fl ipi term 



if I (j)2L term 



1 /" 1 /" d k 

- / ipi{ujcotlibu} +Ki)(fi = - / -—^(t)i{-k) Mi{k) 4>i{k) 

Mi{k) = (ujkCothbujk + Kj) (73) 



/(i k" cIt) -* 

——3 — (pii-k) Ni2L{k,p) <p2L{k,p) 

2p 

Ni2L{k,p) = ^ — ——^.{ujKCsdibujk) (74) 

P - \k2 + ley 



(fl ^ term 



d^k dq 

(27r)3 V 



ipi{-ujcsch.buj)^= I -;-—^—(l)ii-k)Ni^{k,q)(,{k, 



2(7 

Nid^,Q) = ^ n r-^i-u)kcschbujk) (75) 

g^ — [k2 — le)'^ 



£,£, term 



Ij^ {uj + u cothbuj) C = ^ j —-^^^ ^i-k,p)M^ik,p,q)C{k,q) 

f°° dk2 4:pq 

■^«^^'^' '^ = 7-00 ^[p^-(fe2-. 6)2] [g2_(,, + ,,)2] K + c^..coth6a;,) (76) 

hL 4>2L term 

1 /" 1 f d'^k dp dq -^ ^ ^ 

2 / 02L(wcoth5a; + KL)02L = 2 / 7^^ (t>2L{-k,p) M2Lik,p,q) (t>2Lik,q) 
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M2L{k,p,q) 



4pq 



_^ 27T [p2-(k2-ier][q^-{k2 + iey] 



{(jjk coth hojj. + kl) (77) 



hn 4>2R term 



2 / <p2B.i^ + i^r)4'2R = 2 / 7^;^ — — hni-k,?) M2R{k,p,q) <l)2R{k,q) 



M2R{k,p,q) 



dko 



Apq 



2vr [p2 _ (A;2 - k)2] [g2 _ (^^ + ie)2] 



(6jfc + kr) (78) 



^ ^2L term 



^ (w coth ftcj) (j)2L 

Q2Li{k,p,q) 



cPk dp dq 

(27r)2 vr vr 
dk2 



i{-k,p) Q2Lii{k,p, q) 4)2L{k, q) 
Apq 



27r [p2 - (/c2 - ie)2] [q^ - {k2 - ief 



(c^jjt cothfowjt) (79) 



i 4>2R term 



/ ^ (w) (/'2i? 

N2R^{k,p,q) 



d'^k dpdq ^, r \ i.t /r n / /r \ 
YTT-yj ^{-k,p) N2B^{k,p,q) <t>2R{k,q) 

7 '2^[p^-ik2-iey][q^-{k2-ie)^] 



^k 



(80) 



Before taking up the p-dependent terms, we wih consider some simphfication of the expres- 
sions so far. We start with Ai^ and use the integral representation 



OJ + LO coth boj = 



dX oj'^ 
vr A2 + tj2 



/(A), /(A) 



(81) 



1 g—2ib{\—ie) 

The contour for the A-integration is to be completed in the lower half-plane. Using this result 
in (76), we can carry out the fe2-integration. Then we do the A-integral to obtain 



Ms^ = 2tt f{-iu}p)ujp6{p, q) + Apq 



r^ ds 

'o ^ 



(jp. _ g2^(^q2 _ g2-^ (s^ — p2)(g2 — p2) 

f{-iUJq)ujg 



+ 



(s^ — q'^){p'^ — q^) 



(82) 



where ojp = ^/k^ + k'^A-p^, etc., and 



f{-iuj)uj 



2oj 



Y^^=2^ (83) 

The integration over (j)i produces the additional term —Ni^{—k){M.i)~^Ni^, which is written 



out as 



Apq 



-Nu{M^)-^Nu = - /— — - 

^^^ ^' ^^ J 2TT\p^-{k2-ief][q^-{k2 + ief]{sinhhojf{ojcoi\ibuj + Ki] 
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d\ J^ 



vr A2 + w 



MA) 



4pg 



h{X) 



2 "vv [^2_(^2-^e)2][g2-(fe2 + ie)2] 

i\ 

[K/(sm6(A — ie))2 + (A — ie) sin6(A — ie) cos6(A — it)\ 

h{-iKs)Ks 



By carrying out the integrations as we did to obtain (82), we find 



Ni^{Mi)-'N^ = -2n hi-iKp)Kp 5{p - q) - Apq 



+ 



h{-iKp)Kp 



+ 



^ 
h{-iK,)K, 



(jp. _ g2'^(^g2 _ g2 



- + 



with 



h{-iK)K 



[s"^ — p'^){q'^ — p^) [s"^ — q^){p'^ 



{smhbKy{KcothbK + kj) 



Combining (82) and (85) we get 



I^ = 27r {Hp + ujp)6{p -q) + (AM^)pg 



(AM^)pg = Apq 



In a similar way, 



ds 



TT 



{Hs + Ws 



+ 



{Hp + ujp) 



+ 



(Hg + UJg) 



{p^ — s'^){q^ — s"^) {s^ — p'^)[q'^ — p^) {s 



'){p' 



Mr = 2tt (wp + kr)6{p -q) + (AMij)p, 
'•'^ ds " 



{AMii)pg = Apq I 
Jo 



u. 



+ 



UJr, 



+ 



UJ„ 



"^ \_{p'^ — s'^){q'^ — s^) {s^ — p'^){q^ — p"^) {s"^ — q'^){p'^ — q^) 
We follow a similar strategy for A4.2L using the integral representation 

d\ u? 



u coth huj = 
This leads to the expression 



TT A2 + (1^2 



/(A), /(A) = -icot6(A-ie) 



M2L = 27r (iOp coth bup + kl)S{p - q) + Apq / 

^0 



^ 



a;, coth 6a;, 



(84) 



(85) 
(86) 



i7) 



(88) 



(89) 



+ 



LOp coth 6a;p 



(p2 — s2j(|g2 _ g2~) ^g2 _ p2'j^q2 _ p2'^ 

(90) 



Wg coth ftlolq 



is 



'){p — q ) 



We also get a term — A''i2l(— fe)(A^i)~^A''i2L from integration over (/>!. Notice that N12L is the 
same as A^i^ with k2 — >• —k2 and an overall minus sign. These do not affect the final expression 
and we get iVi2L(-^)(-^i)"^^^i2L = Ni^{Miy^Ni^. Combining this with (90), we get 



Ml = 27r {Hp + kl) 6{p - q) + (AM^) 



(AML)pg = Apq 



ds 



vr 



H.. 



pq 



+ 



Hn 



+ 



H, 



{p^ — s'^){q^ — s"^) {s"^ — p^){q^ — p^) (s2 — g2)(p2 _ g2-) 



The expression for Q2L^ simplifies as 



Q2L£, = -4pq 



ds 

TT 



o;, coth 6a;, 



oOp coth 6a;p 



+ 



OJq coth bujq 



(p2 — s2~)^g,2 _ g2^ ^g2 _ p2^^q2 _ p2^ ^g 



'W 



(91) 



(92) 
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which is just the second part of the expression for Ai2L- The integration over cpi also produces 
a new mixing term between (j)2L and ^ given by N-ft{A4i)~^Ni2L which can be simphfied as 



Nl{Mi)-^N: 



12L 



-Apq 



ds 



vr 



h{-iKs)Ks 



+ 



h{-iKp)Kp 



hi-iKg)Kg 



+ 



(s^ — q^){p'^ — q^) 



This expression can be combined with (92) to get 



N2Li, = {Q2Li - Nl^{Ml)-'Ni2L)pg 



-4pq 



ds 



/o ^ 
-(AMl)p, 



H, 



+ 



Hn 



+ 



Ha 



{p^ — s'^){q^ — s^) [s^ — p'^){q^ — p"^) {s^ — <P){p'^ — <P) 



The expression for N2B£^ from (80) can be simphfied as foUows. 

'■°^ dk2 Apq 



N2RC = -4:pq / 



dX a;2 



where 



= 4pq 
= Apq 

= 4pq 

I{k,p) 



_^ 27r [{k2 - ief - p2] [{k2 - ief - q^\ J_^ n X^ + co^ 
dX A2 



27r ^A2 + P (A2 + 6^2) (A2 + ^2) 

dX 1 



27r ^^2 + P 



ijjt., 



I^n 



(A2+a;2) (A2 + 



^7, 



p2 _ q2 



a;2X(fc,p)-a;2x(fc,, 
p2 _ q2 



"^ dX 

TT 



1 1 . f UJp+p 



log 



VA2 + fc2 (A2+w2) 2Tr pup \ujp-p 



(93) 



(94) 



By rewriting l/\/ X"^ + A;2 as we did for the others and simplifying, we can also see that 



N. 



2RS, 



-{AMr) 



pq 



(95) 



(96) 



(97) 



We now turn to the p-dependent terms. 

p p term 

The kernel Mp is cjcoth^w + oj; it is easier to simplify the p p-term after the integration 
over </)i. This yields the kernel Mp = Aip — (cjcsch6a;)2(A^x)^-'^ = H + io. This gives directly 



Mp = 4|/cp 



ds Hs + ojs Co 

2^ wf + ^ 



(98) 



p ^ term 
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This term has cj + w coth buj — Ni^{A4i) ^Nip = H + u as the kernel, so we get 

dk2 2p 2f 



N. 



Pi 



27r p-^ - {k2 - tey ujt 



d\ dko 



4:\k\p 



TT 2tt [p2_(fc2-ie)2](A2+a;2) 
dX A\k\p 



/(A) + h{X) 



°^ ds r"^ dX 4\k\p 

TT 



A\k\p 



vr (A2 + ^2)(A2 + ^2) 

ds iHs +u}s Hp + u}p 



f{X) + h{X) 

fiX) + h{X) 



IT 



UJt 



OJt, 



1 



p2 _ g2 



(99) 



^2L P term 

The kernel for this term is wcoth&cij + C2l — Nip{A4i)^^ Ni2l = H + C2l and gives 

dk2 2p 2\k\ 



N, 



2Lp 



2tt p^ — (/c2 — ieY LtP' 



4|A:|j9 



ds f H^ Hr 



[H + C2L) 



1 



TT V ^2 uD p"^ - S^ 



+ 



C2L 



(100) 



For the first term on the right hand side, we simplified as we did for N^^^ and for the second 
term, we did the k2 integral completing the contour in the lower half-plane. 

(t)2R p term 



This is similar to the <^2LP term, except that we have uj + C2r instead of i:f + C2l- We get 
dk2 ^\k\p 



N, 



2Rp 



2tt [p2 - {k2 - ie)2] OJ 
dXdk2 ^k\p 



^ ,yUJ + C2R) 



CO 



TT 2tt [p2 - (/c2 - ie)2] a;2 ^2 _^ a^ 



C2R 



dko 



A\k\p 



27r [p2 - {k2 - ie)2] w' 



-4|/c|p 



-4|fc|p 



dX 



1 



27r v^A2 + P (A2 + ujD 
2|fc|w2 



2p 

C2i?^ 



cl;^ 



:ioi) 



Finally, we note that in the calculation of the Casimir energy by expansion of the Tr log, 
the first order correction involves the diagonal elements of the type App. In some cases (such 
as AMi, AM^, etc.) Apg is of the form 



Apg = Apq 



ds 



IT 



w.. 



+ 



w„ 



+ 



w„ 



(j)2 _ g2'^(^q2 _ g2^ (^g2 _ p2^(^q2 _ p2^ (^g2 _ q2'^(jp, _ q2'^ 



(102) 
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The p = q limit can be extracted as 

f, r W7_ - W„ /8W\ 1 1 

(103) 



A . 9 I (is 

' TT 



Ws-Wp /'dW\ 1 



2 _ j,2j2 \^ Qp2 J ^g2 _ p2 



We will need this expression with the appropriate W^s to estimate the contribution to the 
energy numerically. 
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